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We investigate the ground-state properties of two-component Bose gases confined in one- 
dimensional harmonic traps in the scheme of density-functional theory. The density-functional 
calculations employ a Bethe-ansatz-based local-density approximation for the correlation energy, 
which accounts for the correlation effect properly in the full physical regime. For the binary Bose 
mixture with spin-independent interaction, the homogeneous reference system is exactly solvable 
by the Bethe-ansatz method. Within the local-density approximation, we determine the density 
distribution of each component and study its evolution from Bose distributions to Fermi-like dis- 
tribution with the increase in interaction. For the binary mixture of Tonks-Girardeau gases with 
a tunable inter-species repulsion, with a generalized Bose-Fermi transformation we show that the 
Bose mixture can be mapped into a two-component Fermi gas, which corresponds to exact soluble 
Yang-Gaudin model for the homogeneous system. Based on the ground-state energy function of the 
Yang-Gaudin model, the ground-state density distributions are calculated for various inter-species 
interactions. It is shown that with the increase in inter-species interaction, the system exhibits 
composite- fermionization crossover. 

PACS numbers: 67.85.-d, 67.60.Bc, 03.75.Mn 



I. INTRODUCTION 



Experimental realization of the atomic mixtures and 
progress in manipulating cold atom systems have opened 
exciting new possibilities to study many-body physics of 
low-dimensional quantum gases beyond the mean-field 
theory. In comparison with the single-component sys- 
tems, mixtures of quantum degenerate atoms may form 
novel quantum many-body systems with richer phase 
structures. Particularly, two-component Bose gases have 
been under intensive studies both experimentally [l], [2, 
HHI1 and theoretically @, 0, H i, Eflj ■ Due to the com- 
petition among the inter-species and intra-species inter- 
actions, many interesting phenomena such as phase sep- 
aration, new quantum states and phase transitions arise 
in two-component cold atomic gases 

SB, 1 1 EH. with 

strong anisotropic magnetic trap or two-dimensional op- 
tical lattice the radial degrees of freedom of cold atoms 
are frozen and the quantum gas is described by an ef- 
fective one-dimensional (ID) model [III EH, [l|[. As 
a textbook example the ID Tonks-Girardeau (TG) gas 
fl4j was observed firstly [H, Il3j . Since then, not only 
the single-component bosonic gas but also ID multi- 
component atomic mixtures have become the current re- 
search focus. In principle, both the intra-component and 
inter-component interactions can be tuned via the mag- 
netic Feshbach resonance, which allows us to study the 
bosonic mixture in the whole interaction regime. Very 
recently, two-species Bose gases with tunable interaction 
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have been successfully produced [3, [j| . 

In general, ID quantum systems exhibit fascinating 
physics significantly different from its three-dimensional 
counterpart because of the enhanced quantum fluctua- 
tions in id mmmHESE^. Thc mcan_fieid the ° r y 

generally works not well for the ID systems except in very 
weakly interacting regime. When interaction is strong 
enough, non-perturbation methods, for instance, Bose- 
Fermi mapping (BFM), Bethe ansatz and Bosonization 
method, have to be exerted to properly characterize the 
features of the system 0, HU, H3, [H, [13, iSl, HEIl^, [H, HI] . 
For a single-component Bose gas, it has been shown that 
the density profiles continuously evolve from Gaussian- 
like distribution of bosons to shell-structured distribu- 
tion of fermions [3(| HH, H2, HH with thc increase in re- 
pulsion strength. Similar behavior of composite fermion- 
ization has been observed in the ID Bose-Bose mixtures 
[27l |28. 29]. For two-component Bose mixture with spin- 
independent repulsions, the normalized density distribu- 
tion of each component displays the same behavior up to 
a normalized constant which is proportional to respec- 
tive atomic number of each component [271] . While the 
extended Bose-Fermi mappings [U [22|, l23j] are restricted 
to the infinitely repulsive limit, some sophisticate meth- 
ods, such as multi-configuration self-consistent Hartrcc 
method [28[ , numerical diagonalization method [2!| , and 
Bethe ansatz method [13, H^, HH, [27j , have been applied 
to study the two-component Bose gases in the whole re- 
pulsive regime. Nevertheless, the above methods either 
suffer the small-size restriction [28l [29| or only limit to 
the integrable models [24), [25|, [26|, [27[ which generally 
do not cover the realistic system confined in an external 
harmonic trap. 
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In this work, we use the basic idea of density func- 
tional theory (DFT) to investigate the ground-state prop- 
erties of the two-component Bose gases trapped in har- 
monic potential in the full interacting regime. It is well 
known that the DFT can handle the system with large 
sizes and will not encounter the exponential wall that 
the above traditional multiparticle wave-function meth- 
ods faced. The effects of the external potential will be ac- 
counted by using the local-density approximation (LDA) 
[H 11, [M, H [33, [H, [H S3, S3, SI . The combination of 
DFT and LDA had been applied to deal with the confined 
ID single-component Bose gas [l9l H3, IH, HI] and Fermi 
gas [33, [H, HE S3, SH, SI successfully. However, the 
two-component Bose gas has rarely been studied. The 
present work will study the ground-state properties of the 
confined Bose mixture and focus on two cases where the 
ground state energy function for the homogeneous system 
can be obtained exactly. In the first case, we consider the 
two-component Bose gas with spin-independent interac- 
tion. In the absence of external potential, this model is 
exactly solvable by the Bethe-ansatz method and thus the 
ground-state energy density function can be extracted 
from the Bethe-ansatz solution [U, [H, HI, H3 • ^ n tne sec " 
ond case, we consider the two-component Bose gas with 
infinitely repulsive intra-component interactions and a 
tunable inter-component interaction. We will show that 
such a model can be mapped into the Yang-Gaudin model 
with generalized Bose-Fermi transformation and thus the 
ground-state energy function can be obtained from the 
Bethe-ansatz solution of the Yang-Gaudin model [43j. In 
the scheme of DFT, the ground-state density profiles are 
obtained by numerically solving the coupled nonlinear 
Schrodinger equations (NLEs) for both the equal-mixing 
Bose mixture and polarized Bose mixture. 

The present paper is organized as follows. Sec- 
tion II investigates two-component Bose gas with spin- 
independent interaction and introduce the method. Sec- 
tion III is devoted to the mixture of two TG gases with 
tunable inter-species repulsion. A summary is given in 
the last section. 



n = 



ti z d 2 



2m,i dx 2 

i=l,2 K 1 
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+ V ext (x) 



(i) 



in which gi (i = 1,2) and g\2 denote the effective 
intra- and inter-species interaction that can be controlled 
experimentally by tuning the corresponding scattering 
lcngthes 01, a 2 and ai 2 , respectively 0, ^t(x) 
(^i(x)) is the field operator creating (annihilating) an i- 
component atom at position x. Here we will consider the 
two-component Bose gas composed of two internal states 
of same species of atoms such that we have the same mass 
for all atoms mi — m 2 = m. When the intra- and inter- 
component interactions are equal {g\ — g 2 = .912 = .9), 
the system is governed by the spin-independent Hamil- 
tonian with the first quantized form: 



JV 



3=1 



d 2 
dx 



2 + Vext(Xj) 



2c^5{a 

3<l 



Xl) 



(2) 



where c = mg/h 2 . The model in the absence of exter- 
nal trap is integrable and can be diagonalized via Bethe 
ansatz method [24 . |26| . 

We first focus on the homogeneous case by setting 
Vext(x) = 0. The effects of the external potential will 
be discussed latter in the scheme of the LDA and DFT. 
Before the application of DFT based on the Bethe-ansatz 
solution, we first give a brief summary of the Bethe- 
ansatz solution for the integrable two-component Bose 
gas. For the eigenstate with total spin S = N/2 — M 
(Q < M < N/2), the Bethe ansatz equations for the inte- 
grable two-component Bose gas take the following form 

sum 



II. BINARY BOSE MIXTURE WITH 
SPIN-INDEPENDENT INTERACTION 



We consider a two-component Bose gas confined in the 
external potential V ex t (x) independent of the species of 
atoms, where the atomic mass of the i-component is rrij. 
The atom numbers in each component are Ni and 7V 2 and 
N = Ni + N2 is the total atom number. The many-body 
Hamiltonian in second quantized form can be formulated 



11 ^ 1 c/2 

1=1 a=l ' 

TV A Z A/ A A 

tan^ 1 t^Z^l = 2nJ a + ^ tan" 1 A ° ~ ^ , (3) 

j=l C < (3^a ° 

where kj and A a are the quasi-momentum and spin ra- 
pidity, respectively, the quantum numbers Ij and J a 
take integer or half-integer values, depending on whether 
N — M is odd or even. In term of the quasi- momentums, 
the energy eigenspectrum is given by E — X)i=i • ^ n 
the thermodynamic limit, i.e., N, L — > 00 with N/L and 
M/L finite, if one introduces density of k roots p(k) and 
density of A roots cr(fc) as well as the corresponding hole 
root densities p h (k) and a (A), the Bethe-ansatz equa- 
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tions thus is simplified to two coupled integral equations 

1 f Q 2cp(k') 
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co-(A) 



-dk 1 



c 2 /4+(fc-A) 2 
cp(fc) 



dA (4) 



2ca(A) 



c 2 + (A - A') : 



■dA'. 



(5) 



The integration limits Q and B are determined by N/ L = 
p = J Q Q p(k)dk and M/L = J B B a{X)d\. 

At zero temperature, the grou nd state corresponds 
to the case with M = [H HE HI, S3, where Ij = 
(N + 1) /2 — j and J Q is an empty set. In other words, 
the ground state corresponds to the configuration cr(A) = 
p h (k) = 0, i.e., no holes in the charge degrees of free- 
dom and no quasiparticles in the spin degrees of free- 
dom. Consequently, the ground states with S = N/2 
are (N + l)-fold degenerate isospin 'ferromagnetic' states 
H, 113, 53, Ej| . Therefore the Bethe ansatz equation for 
the ground state reduces to @, H3] 



p(k) 



9{x) 



1 

27 



1 

27 



2cp(k') 



_ Q c 2 + (k- fc') 2 



2Xg(x') 



2vr 2tt J_ x A 2 + (x - x') 2 



dk'. 



dx' . 



(6) 



where we have made the replacement of variables k = Qx, 
c = QX and g(x) — p(Qx) according to [Hj]. The ground- 
state energy density can be expressed as 

e(p) = |Ve(7) (7) 
where p = £^ . pi is the total density and 

e (7) = 1 1 5(^)^ 2 ^ (8) 

with the dimensionless parameter 7 = f • In terms of 

new variables, the equation p = j_g p{k)dk is rewritten 

as 



A = 7 



l(x)dx. 



(9) 



The expression of ground state energy is identical to 
that of Lieb-Liniger Bose gas [l5| except that now p 
should be replaced by the total density. Here e (7) can 
be obtained by numerically solving the integral equations 
©, © and ® dEEl and we find that it can be fitted 
very accurately by the rational function 



e( 7 ) 



7(l+Pi7+p 2 7 27r2 /3) 
1 + ?i7 + <?27 2 + P27 3 



(10) 




100 



FIG. 1: (color online) Fit function with four parameters and 
two parameters. 



with Pl = 2.05737, p 2 = 0.0688097, q x = 3.16861 and 
q-2 = 0.898297. In Fig. 1 we compare the fit function 
with the exact numerical result, where the exact match 
between them is displayed. In the weakly interacting TF 
regime (7 <C 1) e (7) is approximated as e (7) = 7 and 
in the strongly interacting TG regime (7 > 1) we have 
e(j) — it 2 / 3 such that in these two limits the energy 
density take the form of 



e{ P ) = {y% t 2P ' 

" 6m ' 



7 < 1 
7 > 1 



Alternatively, e (7) can also be formulated by the rational 
functional with only two parameters p = —5.0489470 and 
q = -20.8604983 



e(7) 



7 (l +2?77r 2 /3) 
1 + q-f + p 7 2 



fll) 



which also match well with the exact numerical result 
according to Fig. 1 although the accuracy is less perfect 
than Eq. ifTUf . 

In order to deduce the bulk properties of non-uniform 
systems, we have to combine the analytical solution of 
the homogeneous system with local density approxima- 
tion. Within the local density approximation, one can 
assume that the system is in local equilibrium at each 
point x in the external trap, with local energy per parti- 
cle provided by Eq. ([7]). Therefore the energy functional 
of an inhomogeneous system can be formulated as 





'h 2 


d r 


j dx 


2m 


Tx^ P 



V ext (x)p + pe (p) 



(12) 



where e(p) takes the form of Eq.([7]) which is the exact 
gro und-state energy density of the homogenous system 
[HHlHEEEIil- The first gradient term represents ad- 
ditional kinetic energy associated with the inhomogeneity 
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of the gas that is not accounted for by the "local" kinetic 
energy included in the homogenous energy (3f| H?| ■ Next 
we shall consider the case with external potential being a 
weak harmonic trap V ex t (x) — ^muj 2 x 2 . For the conve- 
nience of latter calculation, we introduce <£>i = ^fp. and 
thus we have /0=|3>i| 2 + | < i ) 2| 2 - In terms of $>i , the energy 
functional can be represented as 



£ = I dx 



Ve, 



$i + Pt (p) 



. (13) 

The ground state corresponds to the minimization of 
the energy functional Eq. (|13p . For the two-component 
bosonic system, the particle number of each component 
is conserved. In order to obtain the ground state from 
a global minimization of £ with the constraints on both 
Ni, we introduce Lagrange multiplier chemical potential 
fXi to conserve N% and p 2 to conserve N 2 . The ground 
state is then determined by a minimization of the free- 
energy functional T — £—p,\ N\ — ^,2^2- As shown in Ref. 
(48|, the minimization of energy functional corresponds 
of the ground state of the modified coupled nonlinear 
Schrodinger equations 

2m dx 2 

h 2 d 2 \ 

— T^+Vext ) $2 + F (p) $2 = P2&2- (14) 

2m dx^ J 

with the normalization condition J dx\$>i(x)\ 2 = Ni, 
where 

F(P) = S- IP* (P)] = [3e (7) - ld,e (7)] . 
op 2m 

In the weakly interacting regime (7 <C 1) and strongly 
interacting regime (7 3> 1) the explicit forms are 



Vext $1 +F(p)$i = 



F( P ) = { 



fJP 



7< 1, 



2 h 2 p 2 /2m 7 > 1, 



while in the middle regime which can be obtained by Eq. 
CD) or Eq. £□}. 

By numerically solving the coupled nonlinear 
Schrodinger equations Eq. (|14p we can obtain the ground 
state density profiles of each component in the full 
interacting regime, which are displayed in Fig. 2 for the 
system with the atomic number Ni — 20 and N2 = 15. 
Here the dimensionless interacting constant U = g/lTno 
and the length x is in unit of I = WU/mw. In order to 
compare the distribution of each component, both of 
them are normalized to one, i.e., rii(x) = pi(x)/Ni where 
J pi(x)dx = Ni. It is shown that in the full physical 
regime the density distribution of two components 
match exactly. That is to say, the density distribution 
of i-componcnt fulfills a simple relation with the total 
distribution 
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FIG. 2: (color online) Normalized ground-state density dis- 
tributions for Ni — 20 (dashed lines) and N2 = 15 (dotted 
lines). From the top to the bottom: U = 0.2,2.0,20,100. 
Unit of rii(x): Ni/l; Unit of a;: I. 



with pt t(x) = p\{x) + p2{x). Similar to the case of single 
species of Bose gas, all of them exhibit the evolution from 
Bose distribution to Fermi-like one with the increase in 
interaction constant. In the weakly interacting regime 
the density profiles display the Gauss-like Bose distri- 
bution, while in the strongly interacting regime Bose 
atoms behave like fermions, which distribute uniformly in 
a more extensive area and the density profiles decrease 
to zero rapidly at the boundary. In order to compare 
the results obtained from the density-functional calcula- 
tions to the exact results obtained by Bose-Fermi map- 
ping method, in Fig. 3 we show the density profiles 
of each component and the total density profile in the 
strongly interacting limit for U = 500. Here each com- 
ponents are normalized to their particle numbers and 
clearly they satisfy the relation of eq. (fl"5"|) . The total 
density profile obtained by BFM method is given by 
Ptg{x) = Yln=a l^nix)] 2 which is also shown in Fig. 3. 
Here 4> n (x) denotes the eigenstate of the single particle 
Hamiltonian of harmonic oscillator. It turns out that the 
results obtained from both methods agree well with each 
other in the full space except that the density profile ob- 
tained by BFM displays oscillations. In the limit of large 
particle number, the differences become imperceptible. 



III. MIXTURE OF TG GASES WITH TUNABLE 
INTER-COMPONENT INTERACTION 

In this section, we focus on the mixture of two TG 
gases with infinitely strong intra-species interaction, i.e., 
<7i = 00 (i = 1, 2) and a tunable inter-species interaction 
.912 =9- In this limit, the original Hamiltonian Eq. (fTJ) 
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FIG. 3: (color online) Ground-state density distributions for 
Ni = 20 (dashed lines) and N2 = 15 (dotted lines) with 
U = 500. Unit of Pi {x): l/l; Unit of x: I. 



is simplified to 
n TG = i <l.r 



i=l,2 



2mj dx 2 
+ 5 *l(x)^(x)* 2 (x)>I<i(x)}, 



+ V ext (x) 



(16) 



with the hard-core constraints ^>l(x)^\(x) = 

= and {^(x), *J(x)| = 1. For the 

mixture of TG gases, it is convenient to introduce the 
generalized Bose-Fermi transformations [49| 



*i(x) 
* 2 (x) 



exp 
exp 
x exp 



nf (z) dz 
nf (z) dz 
(0) dz 



*f» (17) 



*f(x) (18) 



where (x) and \&f (x) are the creation and annihi- 
lation operator at location x for i-component fermions 
and nf{x) — tyf* (x)^f (x). With the above transfor- 
mations, the Hamiltonian (|16[) is mapped into the two- 
component fermionic gas with the Hamiltonian given by 



J i=l,2 



d 2 



2rrii dx 2 
+g*?(x)*P(x)*((x)i>[(x)} 



+ V ex t(x) 



(19) 



The first term on the right side of eq. (fl8|) is intro- 
duced to enforce the fermionic annihilation operators 
&x(x) and &2( x ) fulfilling the anti-commutation rela- 
tions {*f (a;),*!" (x)} = 0. 



Since the bosonic model of (|16p is related to the 
fermionic model of (|19[) by an unitary transformation, 
they have the same energy spectrum. Furthermore the 
ground-state properties which are independent of the sta- 
tistical properties of bosonic system, for example, the 
ground-state density distribution, can be obtained by 
consulting for the corresponding fermionic model. Obvi- 
ously the homogeneous system of (TK)1) with V ext (x) = 
is the Yang-Gaudin model which can be solved exactly 
by means of the Bethe ansatz method (43[. The Bcthe- 
ansatz equations take the form of 



M 



kjL = 2ttIj — 2 tan 

Q = l 



-1 c 



A A —k- 

tan — — 

c/2 



TrJ a 



M 

£ 

(3=1 



tan 



c/2 



(20) 



for the quasi- momentum {kj} and rapidity {A Q } with 
the quantum number Ij and J a . While in the thermo- 
dynamic limit the equations (|20| reduce to two coupled 
integral equations 



p{k) = hj_ 

a (A) = 



dk- 



c'a (A) 



1 

b c' 2 + (fc-A) 2 ' 
I f\ A ,_ -(A') 



-b c 2 + (A - A')' 
Q dk ^ , 



c' 2 + (A - fc) z 



(21) 



with c' = c/2, where Q and B are determined by iV/L = 
f_qp(k)dk and M/L = J B B a(X)dX. If we make the 
replacement of variables x = k/Q, y — A/5, <? c (x) = 
P (xQ) = p (fc) and g s (y) = cr (yS) = a (A) , the above 
integral equations ([2Tj) can be formulated as 



9c (x) 

&> (y) 



1 


1 


27 + 


2tt\ 




£■ 


2ttA c 




1 




7rA s 


L 



dy 



9s (y) 



dx- 



l/4+(x/X c -y/X s y 
9c (x) 



dy' 



1/4 + (y/X s - x/XcY 
9s W) 



1 + (y - y'Y /X 2 S 



with A c = 7 J_ x dxg c (x) and A s = f_ x dyg s (y). Af- 
ter solving the integral equations the energy density of 
ground state can be evaluated by 



£gs 



h 2 N 2 ^ 
2mL 2 Xl 



x g c (x) dx. 



For the system with total density p = pi + p 2 and po- 
larization £ = (pi — P2) / p, the energy density of ground 
state can be parameterized as 

£GS (p ' c) = t( 1/3 +c ^ +f (x ' c) ] (22) 
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with 



f(x,0 



[e( X )-l/3]{l + a(x)C 2 +/3(x) C 4 
-[l + a(x)+/3(x)R 6 }, (23) 



where y = 27/71-, and e (y), a (y) and (3 (y) are given in 
Ref. Ill 51 



P(x) 



-X 



c P X- 
a a X- 



b a 



X 2 + CaX - b a 
a/3X 
X 2 + bpx + c/3 



(24) 



f 7ra p /4, c p = 
1.68894, b a = 
= -1.51457, bp = 2.59864 



with a p = 5.780126, b p = -(8/9)Zn2 
(8/ir)ln2 + 3a p /4, d p = 36 , a a = - 
-8.0155, c Q = 2.74347, 
and C(3 = 6.58046. 

Thus for the mixture of two-component TG gas, the 
energy density e(p) in energy functional (|13[) should be re- 
placed by e GS (p, C)|p-> P (a:),C-»C(aO [H,l39|. With the same 
procedure as the former section, the coupled Schrodinger 
equations take the same formula as Eq. (TTJ} and F (p) in 
the equation about <I>i (i = 1, 2) should be replaced with 

F< (P, = ^-^~P" [(1/3 + C 2 ) + / (x, 0] • (25) 
8m a pi L ' J 

For the unpolarized case, it is reduced to 



F{ P ) 



2*2 



8m 



3p 2 e (x) + P' 



de 
dp 



(26) 



By numerically solving the coupled NLEs, we can ob- 
tain the ground-state density distributions for different 
interacting constants U (U — g/lfuv). 

We first consider the unpolarized Bose mixture with 
Ni = N2, for which we always have p\{x) — p2(x). In 
Fig. 4, we show the density profiles pi(x) for various 
inter-species interactions, which are normalized to the 
atom number N±. When the inter-species interaction is 
zero, no correlation exists between two TG gases and 
pi(x) is given by the density distribution of the single- 
component TG gas, which is identical to the density 
distribution of a polarized Fermi gas with Ni fermions 
fl4j . When the inter-species interaction is weak, the den- 
sity distribution of each component does not deviate far 
from the distribution of the single-component TG gas. 
With the increase in inter-species interaction, the den- 
sity profiles become broader and broader, and continu- 
ously evolve to the limit of strong inter-species interac- 
tion. As shown in Fig. 4, the density profile for U = 500 
already agrees very well with the exact profile in the limit 
of U = 00 except the oscillation peaks. In this limit, the 
density profile for each component fulfils the relation (fT5l) 
and displays N peaks, which can be exactly calculated 
by using the Bose-Fermi mapping method [2l| . 




FIG. 4: (color online) Ground state density distribution for 
unpolarized two-component TG gases with JVi = N2 — 20. 
Unit of pi(x): l/l; Unit of x: I. 
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FIG. 5: (color online) Ground state density distribution for 
polarized two-component TG gases with N\ — 15 and N2 = 
10. Unit of pi(x): l/l; Unit of x: I. 



For the imbalance trapped mixture of TG gases, the 
distributions of each component will display different be- 
haviors in the full interacting regime except in the limit 
of strong interspecies interaction, where the distributions 
display same behaviors only with different normalization 
condition, i.e., n\(x) = ri2(x) with rii(x) — pi(x)/Ni. 
The density profiles of each components and the total 
density profiles for the polarized mixture of TG gases are 
exhibited in Fig. 5a, Fig. 5b and Fig. 5c, respectively. 
The profiles for the case of vanished inter-species interac- 
tion and of infinitely strong inter-species interaction are 
also plotted (oscillating lines) in the figure, which are ob- 
tained by means of BFM method. It is shown that in the 
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regime close to these two limits the distributions agree 
well with the results of BFM method. It is feasible even in 
the situation of finite weak or strong inter-species interac- 
tion. With the increase in inter-species interaction both 
components distribute in more extensive regime and the 
density profiles of i-component evolve from the behavior 
of single-component TG gas with pi(x) — PTc(x)\N = N i 
to the behavior of the fermionized Bose mixture with 
Pi( x ) = wPtg(x)\n=n 1 +n 2 - In Fig. 5d, we show the 
spin density distribution defined as the density differ- 
ence between two components. For U = 0, we have 
pl(x) - p2 (x) = J2n=N 2 \<Pn(x)\ 2 . In the weakly and 
intermediately interacting regimes, the spin density dis- 
tributions display two explicit peaks at the location far 
away from the center of trap and at the locations near 
the middle of trap spin density profiles are almost flat. 
That is to say, two weakly interacting TG gases form the 
shell structure of the coexisted regime of two components 
surrounded by the fully polarized majority component. 
With the increase in interspecies interaction, the peaks 
become smaller and smaller and disappear finally. In the 
limit of strong inter-species interaction, two components 
coexist in the full spaces and the spin density profile ex- 
hibits similar behavior to the total density profile. This 
could be understood because in the infinite-t/ limit one 
has pi{x) - p 2 {x) = Nl ^ N2 ptot(x) according to (TT5j) . 

IV. SUMMARY 

In summary, the ground-state properties of two- 
component Bose gases trapped in ID harmonic traps 



have been studied in the scheme of DFT. On the ba- 
sis of Bethe-ansatz solution, we calculate the ground- 
state density distributions of the confined systems for 
each component by combining the local density approx- 
imation with the exact solution of homogeneous system. 
For the case with spin-independent interaction, the dis- 
tribution of each component is iV, /N of the total density 
distribution in the full repulsive regime. With the in- 
crease in the interaction the density distributions show 
evolution from Gauss-like Bose distributions to Fermi- 
like distributions. For the mixture of two TG gases with 
tunable inter-species interaction, we show that it can be 
mapped into a two-component Fermi gas by a generalized 
Bose-Fermi transformation. Taking advantage of the ex- 
act ground-state energy of Yang-Gaudin model for the 
homogeneous system, we calculate the density profiles of 
the confined systems. For both balanced Bose mixture 
and polarized Bose mixture, the density profiles of each 
component display the continuum crossover from the dis- 
tribution of single-component TG gas with Ni atoms to 
that of composite-fermionization mixture as inter-species 
interaction increases. In the limit of strong repulsion, the 
results obtained from density-functional calculation are 
compared with the results from the Bose-Fermi mapping 
method and it turns out that they match very well in the 
full coordinate space. 
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